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Abstract 

We have observed the south-east (SE) hmb of the Cygnus Loop with Suzaku. Our 
spatially-resolved spectroscopic study shows that a one-ZcTg non-equilibrium ionization 
model represents our spectra fairly well. We find that the metal abundances obtained 
are all depleted relative to the solar values with a positional dependency along the 
radial direction of the Cygnus Loop. The abundances in the very edge of the limb 
shows about half the solar value, whereas other regions inside the Loop show about 
0.2 times the solar value which has been believed as a typical value for the Cygnus 
Loop limb. The "enhanced" abundance in the very edge in the SE limb is quite 
similar to that found in the north-east (NE) limb of the Loop, and thus this is another 
evidence of abundance inhomogeneity in the limb regions of the Loop. The radio map 
shows a quite different feature: the NE limb is in the radio bright region while the 
SE limb shows almost no radio. Therefore, the metal abundance variation in the SE 
limb can not attribute to the non-thermal emission. The abundance inhomogeneity 
as well as the metal depletion down to 0.2 times the solar value still remain an open 
question. 

Key words: ISM: abundances - ISM: individual (Cygnus Loop) - ISM: supernova 
remnants - X-rays: ISM 
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1. Introduction 



The Cygnus Loop is one of the brightest sources below IkeV. It is a nearby (540 pc: 
Blair et al. 2005) proto-typical middle-aged (~10,000 yr) supernova remnant (SNR) located 
near the galactic plane: (/,fe) = (74°, —8°. 5). The foreground neutral hydrogen column density, 
A^H, is estimated to be ~ 4 x 10^°cm~2 (Inoue et al. 1979; Kahn et al. 1980). The low fore- 
ground absorbing material as well as the large apparent size (2°. 5 x 3°. 5: Levenson et al. 1997; 
Aschenbach & Leahy 1999) and high surface brightness enable us to perform a detailed study 
of the soft X-ray emission from the Cygnus Loop. 

The supernova (SN) explosion of the Cygnus Loop is generally considered to have oc- 
curred in a preexisting cavity (firstly noted by McCray & Snow 1979). The X-ray boundary in 
the northeastern (NE) limb of the Cygnus Loop is associated with Balmer-dominated filaments 
which mark current locations of the blast wave (Chevalier & Raymond 1978). Hester et al. 
(1994) studied a Balmer-dominated filament in detail and reported that it was recently (in 
past 1000 years) decelerated from ~400 to ~180kms~^. This rapid deceleration of the shock 
velocity was considered to be a result of the blast wave hitting the wall of a cavity surrounding 
the supernova precursor. Levenson et al. (1999) observed the whole remnant by ROSAT PSPC 
and generated the softness map. They found a soft spatially thin (< 5') shell around almost the 
entire limb of the Loop. Based on the observation, they concluded that the soft shell occurred 
where the cavity walls of the Cygnus Loop decelerated the blast wave. 

ASCA observation of the NE limb of the Cygnus Loop revealed the non-equilibrium ion- 
ization (NEI) plasma condition as well as the depleted metal abundances relative to the solar 
values by a factor of ~5 (Miyata et al. 1994, Miyata et al. 1998). The low metal abundances 
lead them to consider that the plasma in the NE-limb of the Cygnus Loop originated from 
the swept-up interstellar medium (ISM) rather than the eject a- material. Recently, Miyata et 
al. (2007) observed the same region by Suzaku (Mitsuda et al. 2007) Science Working Group 
(SWG) time that enabled us to observe emission lines from highly ionized C and N for the first 
time and obtained depleted abundances: C~0.1, N~0.03, O~0.1, Ne~0.1, Mg~0.2, Si~0.4, 
Fe~0.2. They confirmed the metal deficiency as well as the NEI condition there. Leahy (2004) 
observed the southwestern (SW) limb of the Cygnus Loop by the Chandra satellite and per- 
formed spatially resolved spectral analysis. In his spectral analysis, he divided the elements into 
three groups: 0-group, Ne-group and Fe-group. He found the collisional ionization equihbrium 
(CIE) condition and deficient metal abundances (mean O abundance is 0.22). Levenson et 
al. (2005) observed the southeast (SE) region where a complicated intensity structure is seen. 
They found the metal abundance to be 0.3-0.6 with collisional ionization equilibrium condition. 
Therefore, the low metal abundances may be a common feature except the south east region 
of the Cygnus Loop. 

Cartledge et al. (2004) measured the interstellar oxygen along 36 sight lines and con- 
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firmed the homogeneity of the 0/H ratio within 800 pc of the Sun. We found that the closest 
sight line to the Cygnus Loop was about 5° away. The oxygen abundance they measured is 
about 0.4 times the solar value (Anders & Grevesse 1989). Wilms et al. (2000) employed 0.6 
of the total interstellar abundances for the gas-phase ISM oxygen abundance, and suggest that 
this depletion may be due to grains. Although the ISM near the Cygnus Loop may be depleted 
by a factor of 2 compared with the solar values, the abundances are still much higher than that 
reported at the limb of the Cygnus Loop. It is difficult to explain such a low abundance of 
oxygen in material originating from the ISM. Therefore, the origin of the low metal abundance 
is not yet answered. 

Recently, Katsuda et al. (2008a) found that the northern outermost region in the NE limb 
of the Cygnus Loop showed enhanced abundances: C~0.7, N~0.7, O~0.4, Ne~0.6, Mg~0.3, 
and Fe~0.3 times the solar values. The abundance enhanced region is about 30' x 3' region with 
relatively weak surface brightness. This is the outermost region of the Loop. The inner part 
of it showed depleted abundance. Katsuda et al. (2008b) also confirmed it by Chandra. Based 
on the study of the plasma structure, they concluded that neither a circumstellar medium, 
fragments of ejecta, nor abundance inhomogeneities of the local interstellar medium around 
the Cygnus Loop can explain the relatively enhanced abundance in the region. The apparent 
abundance in the metal enhanced region of the NE limb of the Cygnus Loop can be understood 
as the average metal abundance of the ISM there. Then they studied the possibility that the 
low metal abundance may be due to the contribution of the non-thermal emission since the NE 
limb shows radio bright region (Uyaniker et al. 2004). The Cygnus Loop is one of the radio 
bright sources while it is not uniform. The SE limb region of the Cygnus Loop shows no radio 
intensity. 

We report here the result of the observation of the SE limb of the Cygnus Loop by 
Suzaku showing a similar structure in abundance to that in the NE limb. 

2. Observations and Data Screening 

The observation was performed on 2008 May 12. Fig. 1 shows the location of it by 
a square superposed on the ROSAT image. We employed revision 2.2 of the cleaned event 
data and excluded the time region where the attitude was unstable. Furthermore, we excluded 
data taken in the low cut-off rigidity < 6 GeVc~^. The net exposure time was 18.9 ks after the 
screening. Since the Cygnus Loop almost extends the entire FOV, we subtracted a blank-sky 
spectrum obtained from the Lockman Hole whose observation date was close to that of the 
Cygnus Loop (both are done in 2008 May). We found no emission above 3keV within the 
statistical uncertainties. 

The XIS has a contamination problem (Koyama et al. 2007) that a contaminant is 
accumulated on the optical blocking filter just in front of the CCD. Although the increase of 
the contaminant almost terminates, the low energy response reduces more at our observation 
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time than that at the SWG time. The Suzaku observation of the NE hmb of the Cygnus Loop 
was done for the first time in November, 2005 while our observation on the SE hmb was carried 
out in May, 2008. Due to the reduction of the response at low energy, we used photons in the 
energy range of 0.3-3.0 keV and 0.5-3.0 keV for XISl (back-illuminated CCD; BI CCD) and 
XISO, 3 (front-illuminated CCD; FI CCD), respectively. Figure 2 shows the Suzaku XISl three 
color image. Red, green and blue represent the energy range of 0.31-0.38 keV band, i.e., C VI 
Ka, 0.38-0.46 keV band, i.e., N VI Ka, 0.60-0.69 keV band, i.e., O VIII Ka, respectively. In 
the data analysis, we grouped the spectra into bins of a minimum of 20 counts so that the 
errors are normally distributed to perform a test. Therefore, the degree of freedom (d.o.f.) 
depends on the statistics of the data. 

3. Spatially Resolved Spectral Analysis 

We divided the entire FOV into 28 square/rectangular cells shown in Fig. 2 such that 
each cell has similar statistics. We excluded the lower right square cell since it is almost out 
of the Loop. We extracted spectra from the cells and performed spectral analysis. In order to 
generate the response matrix file (RMF) and the ancillary response file (ARE), we employed 
xisrmfgen and xissimarfgen (version March 2008), respectively. The low energy efficiency 
of the XIS shows degradation caused by the contaminants, however it was taken into account 
in the generation of the ARE file. The latest version of the response includes all the available 
information after the launch (cf. Katsuda et al., 2008a). 

First of all, we applied an absorbed one-ZcTe component of non-equilibrium ionization 
(NEI) model for all the spectra (the Wabs and VNEI model in XSPEC v 12.4.0; Morrison & 
McCammon 1983, Hamilton et al. 1983, Borkowski et al. 1994, Borkowski et al. 2001, Liedahl 
et al. 1995). We noticed that many emission line features are seen while that emissions around 
2 keV or higher are very weak. Therefore, free parameters employed are interstellar absorption 
feature, A^h; electron temperature, kTe] the ionization time, r; the emission measure, EM 
(EM= / ricnudl, where Ue and riH are the number densities of electrons and protons, respectively 
and dl is the plasma depth); abundances of C, N, O, Ne, Mg, Si, S, Fe, and Ni. We set the 
abundance of Ni equal to that of Fe and those of Si/S equal to that of O. This model gave us 
fairly good fits for all the spectra (reduced < 1-3). Maps of the best-fit parameter values are 
presented in Fig. 3. 

Our data clearly show the NEI condition rather than the coUisional ionization equilibrium 
(CIE) condition. If we set logr= 12 where the CIE condition is achieved, we found that 
the reduced increased by ~0.4. Furthermore, we tried to fit data by using a CIE model 
(VMEKAL in XSPEC) which also shows the increase of the reduced by ~0.3 than that of 
the best fit NEI model. Therefore, the spectrum in our FOV surely shows the NEI condition. 

We noticed that the limb regions show higher abundance while the temperature almost 
stays constant. Furthermore, they show radial structure rather than the azimuthal structure. 
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Therefore, we have re-divided the FOV into 10 annular regions from Annl to AnnlO where 
the Annl is the outer most region as shown in Fig. 4. In this way, we performed the detailed 
spectral analysis employing the same model described before. Fig. 5 shows example spectra 
and their best fit models. The reduced values are 1.2-1.7 with 250-440 d.o.f. The detailed 
parameters are summarized in Table 1. 

Fig. 6 shows the variations of various parameters as a function of angular distance from 
the shock front. The value of A;Te is almost constant with a slight increase towards inside. 
The value of r shows an increase towards inside at the very limb region that can be seen in 
the forward shock front (Miyata et al. 1994). We noticed that the metal abundances in the 
outer annuli are much higher than those in the inner annuli. We divide the region into two 
groups: one is the high abundance region [region-H: Annl + Ann2 + Ann3) and the other is 
the low abundance region {region-L: Ann? + Ann8 + Ann9 + AnnlO). The spectra of these 
two regions are shown in Fig. 7 where we can easily see the difference in the line features. We 
also performed the spectral analysis on them by using the VNEI model and summarized the 
results in Table 2. The abundance of O in the region-H is similar to that expected in the ISM 
near the Cygnus Loop (O is about 0.4 solar) while that in region-L is depleted by a factor of 
two and is similar to the average abundance in the shell of the Cygnus Loop (O is about 0.2 
solar). The abundances of C, N and Ne are also high/low in the region-H/L while those of Mg 
and Fe are similar in two regions. The width of the high abundance region is about 3' that is 
larger than that of the PSF of the Suzaku XRT (Serlemitsos et al. 2007). The width measured 
is quite similar to those observed in the NE limb of the Cygnus Loop (Katsuda et al. 2008a, 
Katsuda et al. 2008b). 

If we assume the spherical symmetry of the Loop, we expect that the plasma in the outer 
most region on the apparent image will completely cover the whole Loop. Therefore, we expect 
that the sight line of the region-L is the superposition of the plasma of the high abundance 
region and that of the inside the Loop. If the low abundance comes from the projection 
effect, we can assume that the plasma in the region-L consists of two components: one is a 
high abundance thermal plasma and the other is a very low abundance plasma, resulting an 
apparent low abundance. The low abundance plasma may come from the non-thermal origin. 
Then we applied a combination model of VNEI and SRCUT. The SRCUT model describes a 
synchrotron radiation from a power law distribution of electrons with an exponential cut-off at 
the roll-off frequency, t'roUofr, (Reynolds & Keohane 1999). If we set the metal abundances to 
those obtained in the region-H, the reduced value increased by 250 (744 d.o.f.). If we leave 
the abundance free, we find that the depleted abundance without the SRCUT model gives us 
the best fit result. In this way, we confirmed that the simple SRCUT model does not explain 
the low abundance. This is consistent with the radio result that the SE limb region in the 
Cygnus Loop shows no radio emission. A more sophisticated model will be needed. 
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4. Discussion 



4.1. Abundance inhomogeneities in the limb regions of the Cygnus Loop 

Past X-ray studies have shown that the heavily depleted metal abundance (O abundance 
is 0.1-0.2) is commonly observed in the limb of the Cygnus Loop (Miyata et al. 1994, Miyata 
et al. 1998, Miyata & Tsunemi 1999, Leahy 2004, Miyata et al. 2007, Tsunemi et al. 2007, 
Katsuda et al. 2008a). Among them, Katsuda et al. (2008a) reported for the first time that 
some area in the very edge of the NE limb show high metal abundance (C, N, O abundances 
are ~0.7, ~0.7 and ~0.4). The high metal abundance in the very edge of the SE limb disclosed 
here is quite similar to that of the NE limb, and thus this is another evidence of abundance 
inhomogeneity in the limb regions of the Cygnus Loop. This second discovery may almost deny 
a possibility that the abundance inhomogeneity in the NE limb is just an irregular exception, 
and suggest that such a feature may be observed in other limb regions under some specific 
conditions. 

We here summarize observational facts about the abundance inhomogeneity regions. 
The apparent temperature is relatively low (0.3-0.4 keV) and almost constant with a slight 
increase towards inside. The value of r clearly shows the non-equilibrium condition of plasma. 
The spatial distribution of r also shows an increase from the forward shock to the inside. It is 
quite similar to those obtained in the NE limb regions where the abundances are high (Katsuda 
et al. 2008a). The value of r is large where the abundance is low (inside and bright region) 
both in the NE limb and in the SE limb. This is qualitatively understood that the regions with 
high EM usually show high density and large values of r. If we assume the spherical structure 
and uniform density of the plasma in the SE limb, we can obtain rig to be 0.4 cm~^. With 
taking into account the value of r, we can expect the elapsed time after the shock heating 
to be 2500 years. This rough estimate surely coincides with the cavity explosion. Levenson 
& Graham (2005) observed the SE region just next to our FOV where its surface brightness 
is much higher than ours. They reported the CIE condition at bright region while they also 
reported the NEI condition at the very limb region showing a similar value of r to ours. 

When we see the metal abundance relative to O, the NE limb and the SE limb look 
quite similar to each other with an exception of Fe. In the NE limb, there is a clear difference 
between the abundance enhanced region and the abundance depleted region (Katsuda et al. 
2008a). They can clearly set the boundary between them based on the relative abundances of 
Fe/0 and Mg/0. Whereas, the boundary between them in the SE limb is not so clear, since 
we can hardly see a sharp variation of relative abundances in the FOV. In this way, there are 
some aspects that are different in two regions. However, the abundances are quite similar to 
each other both in the high abundance region and in the low abundance region. Therefore, we 
expect that the origin of the plasma in the NE limb is quite similar to that of the SE limb. 

The low metal abundances may be due to the combination of a thin thermal plasma 
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with normal abundance and a non-thermal plasma. However, our data clearly supported the 
low abundance plasma rather than a contribution of a simple non-thermal plasma. 

4-2. Comparison with the radio morphology 

The Suzaku observation of the NE limb consists of four pointings (NEl to NE4 from 
the south to the north) (Katsuda et al., 2008a). The high abundance region is the very limb 
in NE3 and NE4. Their results indicate that the abundance in the NE4 limb is higher than 
that in the NE3 limb. According to the radio map at 1420 MHz (Uyaniker et al. 2004), the 
brightness temperature, T^, at NEl and NE2 is one of the highest (2-10 K) in the Cygnus 
Loop while that in NE3 is around IK and that in NE4 is less than 0.3 K. Therefore, we find 
that the region where T^ is low can show high abundance in X-ray spectrum. 

Any line of sight through the remnant almost certainly intersects various emitting re- 
gions. However, the line of sight near the limb can intersect the very limited plasma condition. 
Therefore, we speculate that the very limb region in NE3 and NE4 contains only the thin hot 
plasma with high abundance where the radio intensity is weak. We can expect that the limb 
region with weak radio intensity shows high abundance that is much higher metal abundance 
than that in its surrounding region. The radio intensity in the Cygnus Loop is quite strong 
with some exceptions. One exception is the NE limb (NE3 and NE4 and its north) and the 
other is the SE limb. In particular, no radio intensity is reported in the SE limb (Uyaniker et 
al. 2004). We found that both limbs contain regions showing relatively high abundance (O is 
about 0.4 solar) while the inner regions show depleted abundance (O is about 0.2). The NE 
limb and the SE limb are similar in abundance structure while they are quite different in radio 
structure. 

5. Conclusion 

We observed the SE limb of the Cygnus Loop with Suzaku where the surface brightness 
is relatively weak and smooth. We found that the limb region (region-H) shows relatively high 
abundance (0~ 0.4) while that in the inner region (region-L) shows low abundance (0~ 0.2). 
The abundance of O in region-H is similar to that of the ISM near the Cygnus Loop while that 
in region-L is similar to that in the Cygnus Loop observed so far. The high metal abundance 
in the SE limb is quite similar to that of the NE limb, and thus this is another evidence of 
abundance inhomogeneity in the limb regions of the Cygnus Loop. 

The radio map shows a quite different feature: the NE limb is in the radio bright region 
while the SE limb shows almost no radio. Based on the radio intensity morphology on the NE 
limb, we noticed that the high abundance in the very limb region (NE3 and NE4) corresponds 
to the weak radio intensity while the low abundance in the very limb region (NEl and NE2) 
corresponds to the high radio intensity. However, there is no radio intensity reported in the SE 
limb we observed. Therefore, the non-thermal emission does not contribute to the low metal 



7 



abundance both in the NE hmb and in the SE hmb that are left as an open question. 
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Fig. 1. ROSATHRl image of tlic entire Cygnus Loop. Tlic Suzaku FOV in the SE limb is shown in green 
square. 




Fig. 2. Merged Suzaku XISl three color image of the SE hmb (Red: 0.31-0.38 keV band, Green: 
0.38-0.46 keV band. Blue: 0.60-0.69 keV band). The effect contamination is corrected. We divided our 
FOV into 26 rectangles as indicated. 
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Fig. 3. Best-fit parameter maps are shown for various metal abundances. 
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Fig. 4. 10 annuli regions are shown which is employed for detailed spectral analysis. 
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Fig. 5. Example spectra obtained in Ann02, Ann04, Ann06, Ann08 are shown with their best fit models. 



Table 1. Spectral-fit parameters for annular regions 



Parameter Annl Ann2 Ann3 Ann4 AnnS Ann6 Ann7 AnnS AniiQ AnnlO 

feTe[kcV] 0.35±0.01 0.34±0.01 0.38±0.01 0.34±0.01 0.34±0.01 0.35±0.01 0.35±0.01 0.38±0.01 0.42±0.01 0.41 dz 0.01 

C( = Si=S) 0.60± 0.05 0.93 ± 0.06 0.89 ± 0.05 0.56 ± 0.04 0.47 ± 0.04 0.34 ± 0.03 0.26 ± 0.03 0.39 ± 0.04 0.40 ± 0.05 0.26 ± 0.03 

N 0.85± 0.06 1.08 ± 0.06 0.78 ± 0.05 0.52 ± 0.03 0.32 ± 0.03 0.24 ± 0.02 0.24 ± 0.02 0.27 ± 0.03 0.40 ± 0.04 0.24 ± 0.03 

O 0.44± 0.02 0.55 ± 0.02 0.43 ± 0.02 0.31 ±0.01 0.21±0.01 0.18±0.01 0.16±0.01 0.19±0.01 0.22±0.01 0.17±0.01 

Ng 0.69± 0.03 0.92 ± 0.03 0.70 ± 0.02 0.54 ± 0.02 0.38 ±0.01 0.34±0.01 0.30±0.01 0.34±0.01 0.46± 0.02 0.33 ± 0.01 

Mg 0.5±0.1 0.47± 0.09 0.22 ± 0.04 0.29 ± 0.04 0.18 ± 0.03 0.20 ± 0.03 0.15 ± 0.03 0.16='=0.03 0.21 ± 0.03 0.18 ±0.03 

Fo( = Ni) 0.52± 0.04 0.39 ± 0.02 0.27 ±0.01 0.26±0.01 0.19±0.01 0.18±0.01 0.17±0.01 0.19±0.01 0.23±0.01 0.19±0.01 

logC-r/cm"^ sGc) . 10.25± 0.02 10.39 ± 0.02 10.52 ± 0.02 10.54 ± 0.02 10.65 ± 0.02 10.64 ± 0.02 10.7 ± 0.02 10.6 ± 0.02 10.6± 0.02 10.72 ±0.02 

EM[xlO^^ cm~^] 3.92± 0.05 4.14 ± 0.04 5.41 ± 0.05 6.81 ±0.06 lO.lzbO.l 11.1±0.1 12.4±0.1 9.3±0.1 7.0±0.1 11.5±0.1 

X^/d.o.f 392/310 316/254 429/302 457/313 390/302 425/298 433/284 444/257 343/232 578/444 



*Si/S are fixed to O and other elements are fixed to those of solar values. 
The values of abundances are multiples of solar value. 
The errors are in the range Ax^ < 2.7 on one parameter. 
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Fig. 6. Variations of various parameters are shown as a funciton of angular distance from the shock front. 




Fig. 7. X-ray spectra for region-H and region-L. 
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Table 2. Spectral-fit parameters for high and low abundance regions 



Parameter 


high region 


low region 


model 




ViNrjl 


kic[k.c\/\ 


0.34±0.01 


A OA 1 A A1 

0.39±0.01 


C 


A ^A 1 A AO 

0.70±0.03 


A O O 1 AAA 

0.38±0.02 


IN 


n Q 1 _i_n HQ 
U.<3i±U.Uo 


A QQ_1_A no 


(J 


U.44±U.U1 


A 1 A 1 /'I A 1 

0.19 ± 0.01 


AT 

JNe 


0.72±0.01 


A O O 1 A A 1 

0.38±0.01 


J-Vlg 


U.zo±U.U4 


n OA -\- n no 
U.z4 zt U.Uz 


Fe(=Ni) 


0.30±0.01 


0.25 ±0.01 


log(T/cm^'^sec) . 


10.42±0.01 


10.77±0.02 


EM[xlOi^ cm-5] 


0.43 ±0.01 


1.03±G.01 


xVd.oi 


1008/650 


1208/733 



*Si/S are fixed to O and other elements are fixed to those of solar values. 
The values of abundances arc multiples of solar value. 
The errors are in the range Ax^ < 2.7 on one parameter. 
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